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Feshbach's projector technique is employed to quantize 
the electromagnetic field in optical resonators with an arbi- 
trary number of escape channels. We find spectrally over- 
lapping resonator modes coupled due to the damping and 
noise inflicted by the external radiation fleld. For wave 
chaotic resonators the mode dynamics is determined by a 
non-Hermitean random matrix. Upon including an ampli- 
fying medium, our dynamics of open-resonator modes may 
serve as a starting point for a quantum theory of random 
lasing. 

PACS numbers: 42.25.Dd, 05.45. Mt, 42.55.-f, 42.50.-p 



Several recent experiments |^ have demonstrated laser 
action in amplifying random media. The experiments 
use highly disordered dielectrics in which light undergoes 
multiple chaotic scattering. The scattering can make 
light stay long enough inside the material for light am- 
plification to become efficient. A random laser is created 
when the amplification rate exceeds the loss rate due to 
escape from the medium. 

Theoretically, much is known |^ about the sub- 
threshold radiation from random lasers but only a few 
results exist in the non-linear lasing regime. Progress in 
this regime has been hampered by the unusual properties 
of the random laser modes. First, the mode amplitudes 
and mode frequencies in random lasers depend on the 
statistical properties of the underlying random medium. 
Random laser modes therefore must be analyzed in a sta- 
tistical fashion, quite in contrast to traditional laser res- 
onators. Moreover, the character of the modes depends 
on the amount of disorder. For strong disorder, local- 
ization of light may set in and give rise to well separated 
modes centered in different regions of space. By contrast, 
weak disorder leads to a poor confinement of light and to 
strongly overlapping modes. 

Standard laser theory only applies to quasi- 

discrete modes and cannot account for lasing in the 
presence of overlapping modes. Various quantization 
schemes have been proposed to replace the quasi- 
discrete modes of standard theory by quasimodes or 
Fox-Li modes of "bad" resonators. Unfortunately, these 
schemes are not well suited for a quantum statistical 
description required for random lasers. Statistics nat- 
urally enters the random-scattering theory pioneered by 
Beenakker but that approach is restricted to linear 
media and cannot describe lasers above the lasing thresh- 
old. So far, to our knowledge, there is no satisfactory 
scheme for the field quantization in random media. 

In the present paper we develop such a quantization 



scheme for optical resonators with overlapping modes. 
The resonator may have an irregular shape or may con- 
tain weak random scatterers to ensure chaotic scattering 
of light inside the cavity. We employ a technique previ- 
ously applied to condensed matter physics, the Feshbach 
projector formalism. Using that method we show that 
the electromagnetic field Hamiltonian of open resonators 
reduces to the well-known system-and-bath Hamilto- 
nian of quantum optics. Chaotic scattering enters that 
Hamiltonian in two ways. First, the frequency spectrum 
of the resonator modes shows the correlations and "level 
repulsion" typical for wave chaos. Second, the resonator 
mode amplitudes are not damped separately but coupled 
by dissipation. Both effects are related to the spectral 
properties of non-Hermitean random matrices and must 
eventually be included in a quantum theory of random 
lasing. 




FIG. 1. Sketch of a chaotic resonator that is connected to 
the external radiation field via a number of openings. 

We start with the general solution of the quantiza- 
tion problem, and then discuss the application to chaotic 
scattering. For the sake of simplicity we consider a two 
dimensional optical resonator and TM fields; the polar- 
ization vector of the electric field defines the z-axis while 
r = (cc, y) labels the position in the plane. The extension 
to three dimensional resonators and fields with arbitrary 
polarization will be given elsewhere ||] . Several openings 
make for a coupling to the external radiation field. The 
total width W of the openings determines the number 
M of open escape channels at frequency lo, M k, 2yV/A 
with A — 2-Kc/uj. The resonator modes near frequency 
Lo are broadened over a frequency range MAix'/27r, much 
greater than their spacing Auj if Af ^ 1. The resonator 
boundary may have an arbitrary irregular shape. For 
simplicity we assume all walls perfectly conducting. The 
source free Maxwell equations reduce to the scalar wave 
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A(r,i) = 
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for the ^-component A{r, t) of the vector potential. An 
exact quantum description of the total system comprising 
the resonator and the external radiation field is obtained 
by the so-called modes-of-the-universe approach |p|,p^. 
One expands both the vector potential A and the canon- 
ical momentum field 11 = A/c^ in terms of the exact 
eigenmodes '0m(Li;,r) of the Hclmholtz equation 



V'm(t^,r) ^ 0. 



(2) 



A(r, t) ^cj duj q{uj, t)ip{uj, r), 



They are labeled by the continuous frequency lj and the 
integer m; the latter specifies the asymptotic conditions 
far away from the resonator. For example, these con- 
ditions could correspond to a scattering problem with 
incoming and outgoing waves. Then ipm{uj,v) represents 
a solution with an incoming wave in channel m and only 
outgoing waves in all other scattering channels. The 
channel index m may correspond to an angular momen- 
tum quantum number (for a resonator coupled to free 
space) or to a transverse mode index (for resonators con- 
nected to external waveguides). It is convenient to com- 
bine the solutions associated with the different channels 
to an M-component vector iI]{u),y). Then the field ex- 
pansions take the form 



(3) 



where the operators q{Lo) and p{lo) form M-component 
row and column vectors, respectively. Canonical com- 
mutation relation for A{y) and n(r) follow by imposing 
the canonical commutation relations [q„i(ijj),p„(a;')] — 
ih5mnS{Lo — Lo'). The Hamiltonian of the problem is given 

by 

H=lldr [c2n(r, t)' + [V x A(r, t)]'] , (4) 

with A ~ Abz and 11 — liez- The Heisenberg equations 
of motion for 11 and V x A are easily seen to reduce to 
the Maxwell equations. 

The modes-of-the-universe approach yields a consis- 
tent quantization scheme, but does not provide any in- 
formation about the resonator itself. In particular, no 
definition is obtained for the resonator modes, and a 
statistical description of these modes cannot be imple- 
mented at this point. However, further progress is pos- 
sible since the Hclmholtz equation with frequency lo is 
equivalent to a single-particle Schrodinger equation with 



energy E = up- jt? . That link between optics and single- 
particle quantum mechanics allows us to compute the 
mode functions V(w,r) and to define resonator modes. 

The calculation is performed using Feshbach's projec- 
tor formalism |jll|,|l2|. The single-particle Hilbert space 
is decomposed into two orthogonal subspaces associated 
with the resonator and the channel region, respectively. 
The quantum Hamiltonian reads 

+ E /^^ [Wx„m\^>){Xm{E)\+h.c] , (5) 

Am 

with the first two terms describing the decoupled res- 
onator and waveguide, while the last term accounts for 
their coupling. We have chosen a basis in which both 
the resonator and the channel Hamiltonian are diagonal. 
We note that the Hamiltonian (^) is an exact represen- 
tation of the eigenvalue problem (^), even in the regime 
of overlapping resonances. 

The Hamiltonian (H) has been extensively used in the 
theory of chaotic scattering . We employ the formu- 
lation appropriate for scattering through cavities. The 
resonator wave functions 0a (r) are nonzero only within 
the resonator, while the channel wave functions xiE^'^) 
"live" only outside. We require that these functions obey 
Dirichlet conditions along the boundary of the total sys- 
tem (solid line in Fig. 1). The boundary condition along 
the surface separating the resonator from the waveguide 
(dashed line in Fig. 1) is arbitrary save that the total 
Hamiltonian be self-adjoint. The coupling amplitudes 
W\m{E) are given by surface integrals involving appro- 
priate resonator and channel wave functions. 

Diagonalization of the Hamiltonian (||) yields the scat- 
tering states iprn{E, r) with incoming wave in channel m 
only. The modes functions are found from the mapping to 
the Hclmholtz equation, ^{lj, r) = [\/2lj/c]^P{E, r). They 
can be expressed as linear combinations of resonator and 
channel modes 

^(a;,r) =^aA(^)0A(r)+ j dJ (i{Lo,Lo')x{^\v), (6) 



with an Af-component coefficient a\{uj) and an AI x 
M coefficient matrix /3([j,w'). Explicit expressions 
for these coefficients are not needed below. Substituting 
Eq. (^) into Eqs. (||), one obtains the field expansions 

A(r,<) = c^QA<^A(r) + c f dto' Q{lo')x{lo' ,r), 

A 

n(r,t) = i^0*(r)PA + - fdco'xHco',r)P{oj'), (7) 

c ^ CJ 

where we have defined the position operators 
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Qx=J dujq{uj)ax{tj), Q{uj') = J doj q{u;)(3{uj,uj'), (8) 
and the momentum operators 

Px^J dLua{{Lo)p{Lu), P{lu') =. J dujf3\u,u')piio). (9) 

The final step of the quantization procedure is to express 
these operators in terms of photon creation and annihila- 
tion operators. For the resonator modes this is achieved 
by the representation 
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where the matrix U with the matrix elements 
Uxx' = / dr0A(r)0A'(r) 



(10) 



(11) 



specifies the spatial overlap of different resonator modes. 
We note that U is unitary and symmetric, and that it 
only couples degenerate modes, Uxx' ~ S{^\ —^\'), as 
modes with different frequencies have zero overlap (the 
modes are solutions of an Hcrmitean eigenvalue prob- 
lem). The representation (^0|) realizes the commutation 
relations and, at the same time, secures Hermiticity for 
the intra-cavity fields, A = , H = H'' . A similar rep- 
resentation is obtained for the channel modes with the 
replacements ax — > b{ui) and Uxx' U{ujuj'). 

Substituting the representation (|l^) into the field ex- 
pansions (|^) and using the unitarity of we find the 
representation of the intra-cavity fields 
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[a\(f'\{r) - a\(t>l{r)]. 



Substitution of the field expansions into Eq. (| 
yields the field Hamiltonian 

H=^hujxaiax + j duj ujh\^{u))h^{uj) 

A m 

+^E dioWxmiuj) ^a{brniuj) + a{bl^{uj) + h.c 



(12) 
finally 



(13) 



where we defined Wxmi^^) = [c/hy/2uJx]Wxm{E) and 
omitted an irrelevant constant on the right hand side. 
The equations ( p^ and (|l^) are the key results of the 
quantization procedure. The field expansions of the open 
resonator reduce precisely to the standard expressions 



known from closed resonators. However, the field dy- 
namics is fundamentally different as shown below. We 
note that the resonator modes are coupled to the exter- 
nal radiation field via both resonant {a'b, a5^) and non- 
resonant {ab, a'^b^) terms. The non-resonant terms can 
be discarded here since we are not interested in over- 
damping (where mode widths would be larger than or at 
least comparable to the optical frequencies). Our case 
of interest, the case of overlapping resonances, is fully 
compatible with the rotating-wave approximation, where 
only the resonant terms are kept. Then, the Hamiltonian 
( p^ ) reduces to the well-known system-and-bath Hamil- 
tonian ||l^ of quantum optics. It has been argued, that 
this Hamiltonian is valid only for good cavities with spec- 
trally well-separated modes. Our derivation shows that 
such pessimism is inappropriate: the system-and-bath 
Hamiltonian does describe the dynamics of overlapping 
modes, provided the broadening of these modes is much 
smaller than their frequency (so that non-resonant terms 
can be neglected). 

We now discuss the field dynamics and address the 
consequences of chaotic scattering. As a first example, 
we establish equations for the mode amplitudes. From 
the Hamiltonian (|l^) we obtain 

dx{t) = -lioxaxit) -7:J2{WW^)xx'ax'it) + Fx{t), (14) 

A' 

where Fx{t) is the noise operator 

/oc 
du;e-'-it-t'>^Y,WxMcj,to), (15) 

and W the coupling matrix with the elements Wxm ■ The 



equations (14) differ drastically from the independent- 



oscillator equations of standard laser theory, in two re- 
spects: First, the mode operators ax are coupled by the 
damping matrix WW^; second, the noise operators Fx 
are correlated, {fIFx') Sxx', as different modes cou- 
ple to the same external channels (the expectation value 
is defined with respect to the channel oscillators at time 
to). 

A limiting case of Eq. (|lj) is the weak damping regime 
where all matrix elements of WW'' are much smaller than 
the resonator mode spacing Alo. This regime can be re- 
alized either by an opening smaller than a wavelength 
or by the insertion in the openings of partially reflect- 
ing mirrors. To leading order in WW'' / Auj only diag- 
onal elements contribute to the damping matrix, and 
Eq. (jlj) reduces to the standard equation of motion for 
non-overlapping modes |^ . 

For the interesting case of wave chaos the internal 
Hamiltonian can be represented by a random matrix from 
the Gaussian orthogonal ensemble of random-matrix the- 
ory. The eigenvalues lox display level repulsion and uni- 
versal statistical properties. From Eq. (p^, the mode 
dynamics of open chaotic resonators is governed by a 
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non-Hermitean random matrix; we thus encounter an 
interesting connection between the spectral properties of 
open chaotic optical resonators and non-Hermitean ran- 
dom matrices . 

As a second appHcation of our field dynamics we now 
study an open resonator above the laser threshold. We fo- 
cus on a single laser-line and compute the laser linewidth. 
The active medium is represented by Af two-level atoms. 
To simplify the calculation, we assume {i) Af ^ 1 and 
spatially uniform gain , (ii) atomic decay rates that much 
exceed the field decay rates, (iii) exact resonance be- 
tween the laser frequency lu and the atomic transition 
frequency, and (iv) laser operation sufhciently far above 
threshold so that the field fluctuations can be obtained 
by linearization. The Heisenberg equations for the field 
mode amplitudes and the atomic polarization and inver- 
sion take the standard form ^ except for the mode cou- 
pling inflicted by damping and noise. We decompose the 
field into its classical steady-state value and the quantum 
fluctuations, a\ = {a\ + Sa\) enp^—iujt). The steady- 
state conditions take the form = H ■ a where a is an 
A'^-component vector comprising the steady-state ampli- 
tudes a\ (the limit TV ^ oo is taken at the end of the 
calculation). The non-Hermitean matrix 

Haa' = (^A - Cu)Sxx, - n{WW^)xx' + Gdxy , (16) 

depends on the laser field intensity I — J2x\'^>^\^ '^i^ 
gain G = {2SMg'^ hi_){l + Ag'^ I h\\-i ^^^^ . Here, S de- 
notes the pumping strength, g is the atom-field coupling, 
and 7_L and 7|| are the decay rates for the atomic polariza- 
tion and population inversion. The equations of motion 
for the quantum fluctuations follow upon linearization 
around the steady-state solution 
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(17) 



The noise operators T , J-'^ incorporate both field noise 
and noise from the atomic reservoirs. The dynamics of 
5a and Sa'^ is coupled by the 2N x 2N matrix 
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dG / a - a 
'dl[a* -a 



a ■ a 



(18) 



which depends explicitly on the steady-state solution a. 
Equations (p7|), ( p^ reduce the computation of the field 
fluctuations to the spectral decomposition of the non- 
Hermitean matrix C. One easily shows that C has a 
zero eigenvalue connected with the well-known process 
of phase diffusion. The corresponding right eigenvector 
has the form (r, — r*), where r oc a is a right eigenvector 
to eigenvalue of the N x N matrix 7i; the existence of 
r and the corresponding left eigenvector I is guaranteed 
by the steady-state equations, = n ■ a. The phase- 
diffusion coefficient and the laser linewidth Slo can now 
be computed along standard lines We solve the 



equations of motion (|T^) and calculate the Fourier trans- 
form of the stationary correlator {Sa^ {t)6a{0)) . Keeping 
only the zero-eigenvalue contribution in the spectral de- 
composition of £, we obtain the linewidth 



Suj =^ KSlost, 



(19) 



which is larger than the fundamental (Schawlow-Townes) 
linewidth Sujst by the Petermann factor |16|,[17[ 



K^{l\l){r\r) 



(20) 



The non-zero eigenvalues of C will generally modify the 
Lorentzian spectrum and the laser lineshape |18|. A de- 
tailed investigation of these modifications, their statistics 
in a random medium, the pertinent photon statistics [ p^ , 
as well as the generalization to multi-mode lasing will be 
published separately. 
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